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Abstract 

I consider a lattice model of a gauge field interacting with matrix- 
valued scalars in D dimensions. The model includes an adjustable 
parameter a, which plays role of the string tension. In the limit a = oo 
the model coincides with Kazakov-Migdal's "induced QCD", where 
Wilson loops obey a zero area law. The limit a = 0, where Wilson 
loops W(C) = 1 independently of the size of the loop, corresponds 
to the Hermitian matrix model. For D = 2 and D = 3 I show that 
the model obeys the same combinatorics as the standard LGT and 
therefore one may expect the area law behavior. In the strong coupling 
expansion such a behavior is demonstrated. 



1 The model. 



In what follows I propose a lattice gauge model interpolating between the 
Hermitian matrix model Q and Kazakov-Migdal's "induced QCD" f2|. The 
model is very close to standard Wilson LGT. For D = 2 and D = 3 I show 
that the model obeys the same combinatorics as the standard LGT and, 
therefore, if we generally expect area law in the latter then one may expect 
it in the former. Within the strong coupling expansion the area law in the 
model is demonstrated below. The reason for the combinatorial similarity of 
our model to the standard LGT is the symmetry. 

Let us start from the Hermitian matrix model [I| in D dimensions 

Z= [l[D® x e tlV ^ J] e Ntr[ ^ ] . (1) 

x <xy> 

After diagonalizing the hermitian N x N matrices <& x by SU(N) matrices 
U x : <& x — > U x § x Ul and defining for each link < xy > of the D-dimensional 
lattice the SU(N) variables U xy = U^.U y , eq.(^) reads 

Z= [l[d<S> x A 2 (<S> x )e tTV ^ f [] DU xy e NtT[ ** u *y*y u *y ] l[5{I,U f ). (2) 

x "* <xy> f 

where $ x 's are now diagonal matrices, Uf is the ordered product of U xy 's 
along the face (fundamental polygon) / [] and A(<& x ) is the Vandermonde 
determinant. The gauge invariant ^-function is equivalent to a non-poly- 
nomial gauge self-interaction: by definition, the 5-function is a sum over 
irreducible representations r, 

8(I,U f )=J2 d ™rXr(Uf), (3) 
r 

where Xr{U) and dim r = Xr(-0 are characters and dimensions of the r's. 
Using the Cauchy formula, eq.(|) can be expressed as 

N 

J2dim rXr (U f ) = - (U f ) kk }~ N = exp (-JVtr log (/ - U f )) . (4) 

r k=l 
1 In the case of the square lattice, / is the fundamental plaquette. 
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Hence, the {/-dependent part of action takes the form 

OO 1 

S = N te[* ai U*»*vVli +N T,Y, -^W 1 - ( 5 ) 

<xy> / q=l 9 

The first term in (|^) coincides with the action of the "induced QCD" 
[§]. The properties of the model p have been investigated intensively 
and it has been observed that the local Z^v symmetry of the Kazakov-Migdal 
action results in the so called "strong confinement" , or "zero area law" (see 
Ref. for details): the Wilson loop W(C) vanishes unless the spanned area 
A(C) equals zero: 

W(C)=0, A(C)^0. (6) 

If the (minimal) area law W(C) ~ exp(— aA(C)) is presumed, then the 
behavior @ corresponds to an infinite string tension, a = oo. 

We now make the simple observation that the second term in (^), arising 
from ^-function (|3p, breaks local and leads to "deconfinement" 

W{C) = 1 , A(C) ? 0, (7) 

independently of the size of the loop C. Indeed, any product of U xy 's along 
a contour C is equal to product of U/s, where faces /'s fill C. Each such 
Uf = I, due to the 5-function and, hence, Ti< X y>eC U xy = I- It is also clear 
if one recalls that the matrix model (|l|) is actually independent of the U's. 

Formally, this corresponds to zero string tension, a = 0. 

It is natural to try to generalize this model to one interpolating between 
the two limiting cases @ and (0). 

Consider the following simple extension: 

tr log (I-U f ) -> tr log (J — aUf) (8) 

where < a < 1, with a = e Thus < a < oo, and a = oo corresponds 
to the Kazakov-Migdal "strong confinement" (^), while cr = corresponds 
to deconfinement (|7]). 
The action now reads 

S = NJ2 tA^U xy %Ul y ] + NY.T. -tr(C//) 9 . (9) 

<xj/> / q=i y 
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This implies that instead of ([3]) for each face /, we have a factor 

J2^dim rXr (U f ) (10) 

r 

where v r is the sum of all components of highest weight of an irreducible 
representation r = {ni, . . . , n^\. 

N 

u r = ^2n k . (11) 
fc=i 

2 Loop averages. 

If we expect an area law in the standard LGT, then it is natural to expect this 
property in the model Naively, we have checked this already in the two 
limits, (j|) and (0) corresponding to a = oo and a = respectively Now I 
demonstrate for D = 2 and D = 3 that the model with arbitrary < a < oo 
obeys the same combinatorics as the standard LGT and that for large a the 
model obeys an area law. 

Remarkably, a crucial role is again played by the invariance of the term 
tr [QxUxyQyUly] with respect to Zn transformations 

U -> Z N U (U -> UZ N ). (12) 

2.1 L> = 2. 

To calculate loop averages and partition functions we need to average the 
product (over all faces) of the quantities flID"D with respect to the Kazakov- 
Migdal action. The invariance of this action and of the Haar measure under 
transformations flT2"| ) immediately gives selection rules for one-link integrals: 

- for free links (only such links contribute to the partition function): 

J DUe Nt ^^ u "y u \ ri (UV 1 )xr 2 (U^V 2 ) ~ 5 UUU2 (modiV) (13) 

- while for links belonging to a Wilson loop: 

J DUe Nt ^ u ' s 'y u \ ri (UV 1 )xr 2 (^V 2 )ti(UV) ~ 8 Vl+1 , V2 (modiV) (14) 

2 In fact, there is also invariance with respect to arbitrary diagonal SU (N)-\eit and 
-right rotations. 
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These selection rules provide all the combinatorics and, in particular, al- 
low us to calculate Wilson loops at strong coupling (large a). The calculation 
in this case reduces to counting powers of a and to selecting the lowest or- 
der terms in the same way as it can be done in the standard Wilson lattice 
SU(N) gauge theory. For example, for a simple loop with the topology of a 
circle the lowest order non-zero term corresponds to the trivial representa- 
tion for each face outside the contour and to the fundamental representation 
inside. Then, the result is that the total lowest power is equal to the number 
of faces A(C) enclosed by loop C. Thus, we obtain the area law, 

W(C) ~ e ~ aA{c) (15) 

up to a factor independent of a. 

It is worthwhile to remark here that for D = 2 the model becomes exactly 
solvable in the limit of infinite mass scalars. Indeed, following , consider a 
quadratic potential of the scalar field: 

trV($ x ) = N m 2 tr(<S>l) (16) 

After gaussian integration over $'s we have the Kazakov-Migdal action of 
"induced QCD", which is proportional to negative powers of m\. Hence, in 
the limit m\ = oo, the gauge field is completely decoupled from the scalars 
and instead of flTB| ) we have an orthogonality condition for characters 

J DUxniUV^Xr^Vi) = K,r * n d ^ 2) (17) 

and instead of flT4T ) we have the Wigner coefficient 

J DU Xn (U)Xr 2 (tf)tr(U) = D^ n , (18) 

which makes the model to be exactly solvable. The subsequent calculations 
follow the approach suggested in ||. (Full combinatorial details, as well as 
expressions for Wilson loops and partition function on the most general Rie- 
mann surface are explicitly given in 0.) In particular, the result for a simple 
Wilson loop flT5D becomes exact for an arbitrary a: W(C) = exp(— aA(C)). 
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2.2 L> = 3. 

In this case, the symmetry ( |T2"D leads to the following selection rules: 

- for free links: 

J DC/e Artr [*^^ c/t ]x ri (C/V 1 )x J . 2 (C/V 2 )x J . 3 (^3)Xr 4 (^V 4 ) ~ 

~ 5 Vl+VZ)U2+Vi (mod iV) (19) 

- and for links belonging to a Wilson loop: 

j D^e^ tr [^^^ t ]x, 1 (^ 1 )x, 2 (C/V 2 )x r3 (^3)Xr 4 (C/V 4 )tr([/y) ~ 

~ ^1+^3+1,^+^4 (modiV) (20) 

The strategy of further calculation of the Wilson loop in the strong cou- 
pling expansion is straightforward and closely follows the D = 2 case, even 
though the combinatorial details are somewhat more complicated. As a re- 
sult, we again obtain the area-like behavior, 

W(C) ~ ]T X A e-^ c \ (21) 

A(C):dA=C 

where A^'s are coefficients growing slower than exp(aA(C)) f] . The remark 
similar to D = 2 case can be done: in the ml = oo limit the expressions (|19|) 
and (^) become exactly the same as in the pure Wilson LGT. 



3 Conclusions. 

We have a lattice model of a gauge field interacting with matrix- valued scalars 
in D dimensions. The model includes an adjustable parameter a, which 
plays role of the string tension. In the limit a = oo the model coincides 
with Kazakov-Migdal's "induced QCD", where Wilson loops obey a zero 
area law. The limit a — 0, where Wilson loops W(C) = 1 independently of 
the size of the loop, corresponds to the Hermitian matrix model. For D = 2 

3 I do not discuss here the problems related to the competition of coefficients A^'s with 
e - aA ( c ) in the vicinity of RG-stable point, etc. My only purpose here is to demonstrate 
the usual LGT behavior for the model (pi). 
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and D = 3 we have shown that the model obeys the same combinatorics as 
the standard LGT and, therefore, one may expect the area law behavior. We 
have demonstrated this property at the strong coupling. In the massive limit, 
ml = oo, where the gauge field is completely decoupled from the scalars the 
model coincides (up to inessential details) with the pure Wilson LGT. 

Thus, the Hermitian matrix model and the "induced QCD" can be con- 
sidered as the limiting cases of one lattice gauge model. 

In order to better understand how the incorporation of the a term changes 
the matrix model, it is fruitful to look at the D = 1 case (the matrix chain). 
In this case the evolution from a = 1 to a = is very similar to evolution 
from the closed to the open matrix chain, which in turn can be described by 
the action 

M-l 

S = N J2 tr [ ®k$k + i] +aNtr[ ^ M ^} (22) 
fe=i 

where M is the number of sites of the chain, $'s are non-diagonalized. How- 
ever, while in (|22]) a = corresponds to the chain being split, it is not the 
case in our model. The similarity between our model and the model (^) is 
the absence of the non-singlet SU(N) states in the a = limit (the "split- 
ting" point). Obviously, the "splitting" point, a = 0, is a singular one. In the 
language of LGT, this point corresponds to restoration of local Zn, resulting 
in the "strong confinement" . 

The unbroken local Zn symmetry makes it possible to apply the Itzyk- 
son-Zuber formula in the Kazakov-Migdal model || to integrate out the 
LTs reducing the problem to a master field, etc. As a price for solvability, 
the model @ obeys an undesirable "strong confinement" . In our model the 
gauge self-interaction, tr log (I—Uf), breaks the local Z^, but now we cannot 
integrate out the IPs exactly. 
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